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Using lattice field theory techniques, we investigate the vacuum structure of the 
field theory corresponding to perfect fluid dynamics in the Lagrangian prescription. 
We find intriguing, but inconclusive evidence, that the vacuum of such a theory is 
non-trivial, casting doubts on whether the gradient expansion can provide a good 
effective field theory for this type of system. The non-trivial vacuum looks like a 
“turbulent” state where some of the entropy is carried by macroscopic degrees of 
freedom. We describe further steps to strengthen or falsify this evidence. 
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I. INTRODUCTION 


A recent topic of very active interest is to rewrite hydrodynamics as an effective field 
theory [l-4|, with the helds representing the Lagrangian coordinates of the fluid’s volume 
elements. This picture allows the use of well-tested effective held theory techniques to 
investigate huids in the vanishing viscosity limit, a limit where the very dehnition of hydro¬ 
dynamics is currently ambiguous js, 6|. An advantage of the held theory approach is that 
thermodynamic concepts like huid isotropy and entropy conservation can be represented as 
symmetries. 

Phenomenologically, such a theory can be applied to a wide variety of settings, from 
superhuid helium to cosmology ^ to quark-gluon plasma Is I. From a theoretical point of 


view, it allows us to access a hitherto unexplored region |9l-lll|: One where the mean free 
path is small enough to neglect all dissipative ehects but where the microscopic number of 
particles is not “large”, so microscopic decorrelation (molecular chaos, or, equivalently the 
large Nc limit in an AdS/CFT jo] setting) does not apply. Hence, fherma/huctuations excite 
hydrodynamic degrees of freedom which subsequently evolve non-linearly: When viscosity is 
so low that “typical” sound waves, of frequency ~ T and amplitude comparable to a thermal 
huctuation, Ap/ (p) ~ Cy/T^ (where p is the energy density, Cy the heat capacity and T the 
temperature), survive for a time much larger than the thermal scale, ~ 1/T, Kubo’s formula 
needs to be renormalized to account for the energy-momentum carried by the sound waves. 
This is a reasonable physical interpretation of applying “quantum” concepts, such as the 
dehnition of observables in terms of functional integrals, to something so quintessentially 
classical as a “perfect huid”, and might be used to demonstrate that the existence of a 
quantum limit on viscosity is plausible from hydrodynamic arguments alone [h|. 

Works such as js, 12, 131 have established a consistent set of techniques of generating 
EFT terms in increasing order of derivatives respecting the fundamental symmetries of huid 
mechanics, with recent dissipative applications [l4, 15| becoming possible. However, it is 
well known that for an EFT expansion to be complete, the right vacuum of the theory has 
to be known and expanded around. The presence of turbulence in classical huid dynamics 
suggests that vacuum choice around the hydrostatic hunt, employed in R Q Q. is not 
necessarily justihed, an assumption also discussed in [2[. 

Figured] illustrates how this could happen: While the microscopic degrees are traced over. 
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they give a contribution to the free energy, dimensionally set at Tq , where Tq is the cutoff 
scale for microscopic degrees of freedom, which needs to be dominant for the hydrostatic 
limit to apply. If the heat capacity at constant pressure of the system is high enough 
w.r.t. the compressibility (note that the planar limit explicitly excludes this, since there the 
heat capacity trivially vanishes) thermal fluctuations will seed waves and vortices which do 
not go away but will interact non-linearly, modifying the hydrostatic vacuum. In terms of 
fundamental statistical mechanics, this dynamics can be explained by the entropy carried 
by the microscopic DoFs (strictly determined through the temperature Tq and the equation 
of state) vs the macroscopic DoFs (determined by the attice spacing and the functional 
integral of the helds). When the latter becomes a sizeable fraction w.r.t. the former, part of 
the microscopic entropy will be shared, under conditions of equilibrium, with macroscopic 
excitations. If this conhguration dehnes a minimum in the free energy, the hydrostatic 
vacuum is unstable and a non-trivial phase diagram results. 

Since the presence of turbulent flow destroys the symmetries of the hydrostatic limit 
(at least compressibility, but in principle homogeneity and even isotropy), the dynamics 
described above, if it exists, must give rise to a phase transition. Figure [T] assumes a 
hrst-order phase transition, where the two minima coexist at a given temperature. Second- 
order phase transitions, critical points and cross-overs are also possibilities, set by the fluid- 
dynamical free energy in the standard way Q. 

The purpose of this work is to test this assumption using Lattice Monte Carlo simula¬ 
tions, the only known way to examine a quantum or statistical theor y in dependently of its 
perturbative structure, continuing the preliminary analysis shown in 181. 


We would like to emphasize, in case of possible misunderstandings, that the calculations 
done here are not real-time simulations of classical solutions such as 19| , but a simulation 


of a static “vacuum” whose partition function is evaluated by a functional integral, ana. 


gous to lattice QCD simulations (and in fact using the same numerical techniques 20, l21|L 


o- 


appropriate to describe, beyond any linearization or perturbative expansion, both quan¬ 
tum/thermal fluctuations and their response. Our results are therefore not an evolution 
of the theory from initial conditions, but a description of the “ground state” of the theory 
which, together with the Lagrangian, dehnes the equations of motion. If this ground state 
is the hydrostatic limit (a hnite temperature static huid), hydrodynamics will arise as an 
effective theory with the gradient dehning the coarse-graining scale. If it is not, just as in 
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FIG. 1. A schematic illustration of a possible vacuum structure of the theory, in terms of a free 
energy including both microscopic degrees of freedom and macroscopic collective excitations. The 
left side is the ’’hydrostatic limit”, where expansions such as 3(] lead to a well-behaved effective 
theory. The right side is a possibly turbulent vacuum where a fraction of the entropy goes into 
microscopic degrees of freedom. In this figure the free energy corresponds to the coexistence phase 
of a first order phase transition, but of course other configurations between the two vacua (second 
order, cross-over, critical points and so on) are equally possible 


Image on the right from SCIDA 


El 


usual QFTs, vacuum effects will introduce additional correlations which will not be visible 
to any order in the gradient expansion (a good analogy here is the role of instantons in QCD 

3 ). 


II. THE THEORY 


A. Lagrangian description of hydrodynamics 


Three-dimensional ideal hydrodynamics with no chemical potential (any “particles” are 
balanced by “antiparticles”, so the net density of each conserved charge vanishes in each 
volume element) can, in the Lagrangian picture, be described [l, 2| in terms of three fields 
0/=i,2,3^ which physically correspond to the x,y,z coordinates of the comoving frame (La- 
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grangian picture) with respect to the lab frame (in which the Eulerian picture is dehned). 
Note that unlike scalar helds used in quantum held theory, cj)^ have dimensions of one unit 
of spacetime. 

The choice of is of course not unique, as a perfect huid is homogeneus, and in 

its comoving frame, invariant under rotations and rescalings. This symmetry restricts the 


Lagrangian to the form 


,a 


L = F{B)=T^F{det{Bij)) , Bjj = . 


( 1 ) 


The function F{B) is left arbitrary, as it corresponds to different equations of state for 
the huid. We note that the symmetry automatically determines the interactions, and the 
coupling constant of the theory, which therefore does not need to be hxed on the lattice. We 
also note that any excitations of the theory can actually be regarded as goldstone modes 
of one of the broken symmetries included in this Lagrangian j^. This makes it imperative 
to distinguish dynamical symmetry breaking ehects from lattice artifacts. The next section 
discusses some ways this was done, but, while we believe the evidence is very good that the 
results presented here are not artifacts, it is not conclusive. For this reason our results are 
labeled as “indications.” 

Dimensional analysis makes it apparent that the F{B) should be dehned in terms of 
an energy scale Tq, which can be identihed with the characteristic scale of the microscopic 
degrees of freedom, or equivalently the “microscopic temperature” of the system. Note that 
this tells us only about the density (and huctuations of it) and is in general different from 
the mean free path of the interacting theory, which in the ideal hydrodynamic limit goes to 
zero. Normally, for the Boltzmann equation (and, in AdS/CFT, classical gravity) to make 
sense, Tq needs to be much smaller than the mean free path. This work takes the opposite 
limit, since the mean free path is zero and Tq is hnite. 

If Eq. ([1]) is used to build a partition function, the “effective Planck’s constant” becomes 
dimensionful, as the “microscopic gradient” entirely factors out of the Lagrangian. Expan¬ 
sion around Tq, therefore, is potentially very different from the gradient expansion since 
“non-perturbative” contributions (~ exp[—l/(a;To)], where x is some distance scale) never 
quite go away. Mathematically, this is analogous to non-Abelian gauge theory, where the 
action ~ gyl^F"^ and likewise one cannot generally expand in powers of qym except in some 
rigorously dehned limits (the connection between non-perturbative physics in Yang-Mills 
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theory and turbulence could in fact be more extended [23|; indeed, it has long been clear, 
and this paper may partially conhrm, that Wilson loops and vortices share deep similarities). 

It is straightforward to show that the classical expectation value of the energy-momentum 
tensor corresponding to the Lagrangian in Eq. ([1]) is that of ideal hydrodynamics 24 1 


{T^^u) = {p + p)u>^u’' + 


( 2 ) 


(using the “mostly plus” metric) and hence this is simply an unusual reparametrization of 
ideal hydrodynamics. The energy density and pressure in this notation are 

dF 


p = -F{B) , p = F{B) - 2B 


dB 


( 3 ) 


Hydrodynamic flow is dehned as being perpendicular to any gradient of the 0^: u^dncj)^ = 0 


This, and = — 1, unambiguously give 


= 


6\/H 


eijKda(i)^ d0(i)'^ d^(f) 


K 


( 4 ) 


We can also show that d^{-\fBu^) = 0. By inspection, without any conserved charges (those 
are examined in j^), one can identify 


= gT!^ 


( 5 ) 


with the microscopic entropy. Using the Gibbs-Duhem relation, then, the temperature will 
be 

p + p 


T = 


= Tr 


^fB{dF/dB) 


( 6 ) 


s 9 

Note the presence of p as a free parameter. This is the microscopic degeneracy, an intrinsic 
property of the system. In the planar limit of Yang Mills theories, it tends to inhnity 
faster than any other constant of the system, making the fluctuations discussed in this 
work irrelevant {Tq/q —)■ 0 and any microscopic fluctuation gets distributed equally, by 
equipartition, between a “large” number of degrees of freedom). 

The likely ubiquity of non-perturbative effects in this theory can be demonstrated by ex¬ 
amining the “vortex” degrees of freedom, as was done in j2|: Naively, a vortex can be treated 
as a non-topological inhnitesimal perturbation, a phonon ttt equivalent to a soundwave vr^. 
However, in a hydrostatic background, vortices do not propagate, yet carry arbitrarily small 
amounts of energy and momentum and interact in ways constrained by the symmetries of 
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the Lagrangian {tiltit -H- til'll and vr^ —)■ til'll are possible) Thus, fluctuation-driven vor¬ 


tices can become stable, and their interactions could overwhelm the vacuum state. In 


25| 


this interaction is demonstrated to have attractive componen 


of “vortex condensates” forming in the vacuum. However, as 25(] treats vortices as source 


^ leading to the likelihood 


terms, the equilibrium state of a liquid subject to fluctuations and its dependence on Tq 
cannot be ascertained through such a perturbative expansion. The alternative analytical 
approach, to deform the theory in the infrared [l, 11 1, is liable to give an incorrect vacuum 
since the deformation breaks the symmetries of the Lagrangian. This leaves the lattice as 
the only avenue to investigate vacuum properties consistently. 

The theory formulated via Eq. [T]can be put on the lattice 3, 21] in the usual way, via 


\nZ = / V(f)iexp\i / d xL + J(j)jj / d^^exp 

^ ^ lattice+Wick'^ 

Throughout this work we use the ideal Bose gas EoS, 

F{B) = . 




( 7 ) 


( 8 ) 


This can be easily generalized to any monotonic EoS (without phase transitions), for example 
an EoS, say, fitting the QCD cross-over [ll|. The ideal gas, however, is a good testing 
laboratory as all of its parameters are very simple 


= . {p) = k> , = 

\Hi 3Q , \F/ 3 , \ / ^3 


( 9 ) 


Again, note that the starting point in this formalism requires the identification of s oc \/H 
and, from the formulas above, the following relation: 

4/3 


(P) = 


TT^ / 45 


d' 


(s) 


4/3 


( 10 ) 


30 \2n^) 

where the constant of proportionality (~ 0.987) can be absorbed into the ratio of scales 
(aTo)^, giving the above F{B). 


B. Lattice implementation 

Since centered differences lead to two disjoint, interwoven lattices, our lattice derivatives 
are approximated by an average of all (eight) one-sided differences per hypercube. All 
















quantities derived from the field derivatives 9^0^ (e.g., -B/j, are therefore situated 

at the centers of the hypercubes. To handle the periodic nature of the toroidal lattice, we 
subtract the hydrostatic background (when 0'^ = x^) and deal with the “shifted” coordinates: 

71^ = (f)^ - -)■ 0 „ 0 ^ = daTT^ + ISl . ( 11 ) 

The TT'^ are the fluid phonon helds. 

Because we expect extended structures (e.g., vortices) to arise, we use HMC updates. 
The required variation of the action with respect to local held values is given by: 

6Vb < 5 ( 0 ^ 0 -^) 

6(j)^{x) 6y/B S{da(j)-^) 5(j)^{x) 


E 




dF 

dy/B 


S^'^6{y — x± /i/2 ± ij 12 ± (j/2) 




B^k 


-fiiyaa 


da(j) 


K 


y-al2 


y+al2 


The updating algorithm and the calculation of observables has been implemented in C code 
(including the ranlux random number generator 2^), with multicore parallelization via 
OpenMP. 


Interpreting the lattice data is complicated by the fact that this is a manifestly non- 
renormalizable theory: The continuum limit is one where Tq —)■ cxo, and the approach to it 
is obviously divergent for any dimensionful quantity, though not for dimensionless numbers. 

This however is physically reasonable, since the regulator in this theory can be interpreted 
as the approach to the physical microscopic scale. Varying aTo, the lattice spacing in units of 
the microscopic scale, we are considering different systems whose ratio of the characteristic 
size of microscopic excitations to collective excitations is different. This, in this context, 
is the only free parameter of the theory. (However, another possible one arises via the 
“macroscopic” temperature, the period in the time direction of the hnite temperature lattice, 
and the influence of this should be investigated as well.) 

The continuum limit would be an inhnite lattice of a given aTo. This limit is obviously 
not reachable, or even approachable by the hnite computing resources used so far in this 
project. Our results thus far, however, show that we may be closing in on such a limit. 
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FIG. 2. The energy density as a function of oTq and the lattice size for ensembles both with and 
without collective “sound modes” present. Left; Average energy density in lattice units. Right: 
ratio of average energy density and the 4/3-power of the lattice-averaged entropy density, (s). 


III. RESULTS 

A. Evidence for coexisting phases 

The overall energy density p and its relation to entropy density s are explored in Fig. |2l 

Beyond a “critical” aTg (~ 0.47 for the lattices), two vacua are clearly visible. The 
upper vacuum’s entropy and energy has a monotonic power law dependence on aTo, and can 
be readily identihed with a state where the bulk of the entropy is carried by microscopic 
degrees of freedom. The lower vacuum has a nearly constant microscopic entropy with 
respect to aTg, with a nearly constant ratio of energy to entropy density, and can be identihed 
with the vacuum where entropy and energy density has gone into collective (“sound”) modes. 
One can see in the plot on the right that corresponding entropy-density correlations have 
managed to maintain an enhancement in (p) / when such modes are present. 

To interpret this data we note that “low” and “high” entropy states refer only to the 
microscopic entropy, since we have no way, at present, of measuring the entropy contained 
in the macroscopic perturbations (see, however, the Discussion section, specihcally around 
Eq. [T7|l . Thus, the enhancement of (p) / is in line with the hypothesis that macroscopic 
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degrees of freedom in the “low entropy” phase carry a significant fraction of entropy, as well 
as energy density. 

The locations of the vacua show a small dependence on the lattice size, indicated by the 
different symbols in the figures, with the jump to the 20*^ and 32^ lattices thus far showing 
little difference from the 16^. 

Transitions between vacua have also been observed, but appear to be exceedingly rare. 
Close to the “critical” point, we have observed single transitions from the lower entropy 
state to the higher one; and well above this point, we see single tranisitons from the higher 
entropy state to the lower one. Hence, thus far we are not able to verify the scaling of these 
transitions with system volume, which, if its rate T went as 

Ti^s ~ exp [-|si - S 2 I H] (13) 

would be a clear indication of hrst-order behavior. 

10° 
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FIG. 3. Plot of the time component of the flow as a function of uTq. 
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Figure [3] compounds these indications by examining the average of the (Euclidean) time 
component of the flow Uq. While spacelike components (m*) are on average zero if rotational 
symmetry is not broken (and by observation it is not), (uq) shows a clear non-zero value 
rising steadily towards the hydrostatic limit (^ 1). Along the way, however, the lower- 
entropy vacuum arises and transitions to it result in substantially increased spatial velocity 
fluctuations, lowering (uq). 
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This flow is not fully explainable in terms of perturbations, as can be seen by the fact 
that, for higher values of oTq, the average value of the flow is multi-valued, something only 
possible in the presence of a first order transition: In this regime, the vacuum where uq 
is lowest (and hence where fluctuations in Ui are higher) has an average expectation value 
of collective excitations, which agrees well with the idea that this vacuum is dominated by 
fluctuation-driven turbulence. 

In the future, correlators of flow observables such as vorticity 

Cp = ^ {p + e)u^dx^ (14) 

or of the tensor shear 


= {Uf,Uy + , (IS) 

(the average values are fixed to zero by symmetry) can be examined in the two regimes to 
see what difference is there between the flow structures of the two vacua. 


B. Evidence for “sound” modes 

Figure m shows the space-time dependence of the (hypercube averaged) phonon correla¬ 
tors. On this ensemble, one can clearly see a nearly time-independent longitudinal mode in 
(7r*(0)7r*(a:)). The absolute magnitude of the fields does not matter, but such correlations 
will clearly enter into the derivatives (and hence etc.). “Transverse” modes (e.g., 
sinusoidal x-dependence in (7r^(0)7r^(a;))) are also sometimes seen, sometimes in combination 
with the longitudinal ones. 

In Fig. [5] we plot corresponding relative entropy-density correlators. Looking at the plot 
on the left, one can see that such correlations fall off into noise rather fast for ensembles 
without sound modes (i.e., from the higher energy branch), whereas for the lower energy 
ensemble, they persist across the lattice. It is straightforward to relate these persistent 
sound modes to the average value of the flow component in Fig. [21 However, as the next 
subsection explains in detail, in a finite system, care must be taken to interpret these degrees 
of freedom as mere sound perturbations. 
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FIG. 4. Phonon correlators of longitudinal “sound” modes. 



t'/a x'/a 


FIG. 5. Relative entropy-density correlators. Left: results for two 16^, aTg = 0.5 ensembles - one 
with a “sound mode” and one without - as a function of time separation (dotted symbols are for 
negative values). Right: the spatial dependence from the sound-mode ensemble. 


1. Interpretation in terms of kinks and domain walls? 


Another feature of first-order phase transitions is the presence of kink solutions in the 
mixed phase: If two global vacua are possible, a translationally dependent interpolation 
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between them is also a minimum. Such objects, in conventional phase-transition physics, 
are called domain walls, solitons, instantons (if time-dependent), calorons (if at hnite tem¬ 
perature), and so on. 

In the context of classical hydrodynamics, a natural implementation of this concept is the 
relativistic static shock: long-lived shocks, while generally unstable against small perturba¬ 
tions, arise as classical solutions to the hydrodynamic equations of motion. Relativistically, 


the Taub adiabat is the most well-known example 


24 


271. l28| . If the equation of state is not 


single-valued, these shocks can be stationary [2^ |28| and, in presence of a surface tension, 
stable. Such topological objects form the surface of objects such as bubbles, which charac¬ 
terize the mixed phase of a hrst-order phase transition. In a sudden cooling, they can also 
remain in the “colder” phase, where symmetries of the hot phase are broken. 

In the context of our hypothesized hrst-order phase transition between the turbulent and 
the hydrostatic vacuum, the kink should appear as a long-lived shock between a hydrostatic 
and a boundary layer. The shock’s width could be interpreted as a “lower limit” to Prandl’s 
boundary layer 29|], which would reinforce our contention that quantum hydrodynamics can 
be seen as a lower limit of viscosity. 

Indeed, Fig. Hlshows a correlator of the how held tt^ (See Eq. [TT]) where a space-like mode, 
stable in time, can be seen. The stationary nature of this mode precludes it from being a 
perturbative expansion, but it is somewhat expected if the dynamics of the vacuum is driven 
by vortices, since vortices are stationary j^. 

Its space dependence could be seen as an indication of a topological structure, but for 
conclusive evidence of this we would need to study three point functions of correlators such 
as Eqs. (HI and ina, to conhrm that the vacuum structure around a kink interpolates 
between the two phases. 


IV. DISCUSSION 

A. Effects of the lattice formulation 

As mentioned earlier, we cannot remove the hnite lattice spacing, but we can alter its 
effect by choosing different discretizations for the lattice derivatives. More extended hnite 
differences are certainly possible and we could even “take a step backward” and move to 
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centered differences to explore the effect of finite a. (Much older results 30| , from the start 
of this project, using centered differences, indicate a “critical” point around aTo ~ 0.7 on 
lattices, close to where we see one now on lattices of the same size.) 

Since the collective modes which we see on our lattices tend to £t with one wavelength 
within the hnite, periodic box (see Fig. Hj), we cannot help but wonder whether these are 
simply symptoms of the toroidal nature of the volumes they occupy. It may be worthwhile 
switching to Dirichlet {tt^ = 0) or Neumann {diir^ = 0) boundary conditions to see how 
they behave there. But we point out that these modes also arise from the particular action 
in question (see Eq. |H]) and that they more readily appear for larger volumes, making us 
believe that they (and their effects) will persist in any “thermodynamic limit” {V —)■ oo). 

Another consequence of using the lattice formulation is that we necessarily have a Eu¬ 
clidean signature (see, e.g., (mq) in Fig. [3]). One might imagine that, close to the hydrostatic 
limit ((mq) ^ 1), we could map our Euclidean hydro to a Minkowski one: e.g. 


= (±Vl - = (±1/7.lO -> = (-TiTv) > 


(16) 


but we find that we are still quite far from this limit when the collective modes (and larger 
velocity fluctuations) appear, driving the limit beyond our control. 

Before moving on we make one more observation about the two “phases” we observe: 
larger spatial volumes drive the higher energy state to even higher energy and entropy and 
the lower energy phase to lower energies, but larger time extents have the opposite effect 
for only the lower-energy (collective-mode) phase. We may be in a position to explore 
the “macroscopic” thermodynamics of these collective excitations. The interaction measure 
(/ = p — 3p) can expressed as 


+ /o = 


Tmacro d 111 Z 


V din a 


T 

-L 't 


macro ^ I d ln Tq 


V 


din a 




, dlnTn , , , 

= -4(1 + ^^) (p). 


.1 , ( 17 ) 

d Ina 

where we used Tmacro = l/(aAi) and V = and Jq represents the zero-temperature 

(divergent) part which must be subtracted off. So we see the interaction measure of collective 
excitations is related directly to the average microscopic energy density. From here, the 
integration method (see, e.g., |3]|) can be used to arrive at the macroscopic pressure, energy 
and entropy density. 
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B. Further developments 


The main barrier to our results not being more conclusive are the computational limits. 
This paper is invariably a series of several ones, with higher statistics runs giving a more 
coherent picture of the situation. 

The availability of more statistics will allow us to understand if the transition between 
the two vacua follows the volume scaling law expected from a hrst-order phase transition 
(i.e., Eq. [13]). We would regard this as definite proof hydrodynamic fluctuations trigger a 
first-order phase transition in the vanishing viscosity limit. This proof is at the moment 
lacking. 

In Section IIII Al we have proposed investigating turbulence using vorticity and flow cor¬ 


relators. Such correlators are well-studied in classical turbulence 


35| . In particular, if our 


vacuum is similar to classical hydrodynamic turbulence, we would expect that they fall as a 
power-law, with an exponent related to turbulent spectral exponents (such as the well-known 
Kolmogorov spectrum) calculated in semiclassical theories. It would be very interesting to 
investigate this further. Studies such as these are instrumental in formulating the effective 
Lagrangian of the colder phase. If our interpretation is correct, turbulence results in the 
symmetries of hydrodynamics being broken. Thus, the effective Lagrangian of the cold phase 
will look like 


= F{B) + J2<^^iTo)f{Bn) + J2UTo)g{Bu) + (18) 

ill 

where coupling constants a*, 6*, q break the conformal, isotropy, and homogeneity conditions 
respectively. For a true dynamical symmetry breaking effect rather than a lattice artifact, 
these terms must arise as a phase transition rather than a smooth cross-over, something 
which we seem to observe. Furthermore, operators of dimension D should be Tq if the 
theory examined here is indeed an effective theory valid at scales much smaller in momentum 
space than Tq. 

This study could lead to a full non-perturbative renormalization group analysis of hy¬ 
drodynamics as a field theory. Naively, this field theory is non-renormalizeable, something 
seen both on dimensional grounds (appearance of operators of dimension higher than four) 
and on physical grounds (the coarse-graining scale in hydrodynamics is a physically ob¬ 
servable quantity). This question is directly connected to the question of whether “ideal 
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hydrodynamics” with the coarse-graining scale going to infinity (in other words the limit 
r]/s 0) is well-defined. In general, for an observable X, iff lim^^o {X) ~ {^o)termostatic’’ 
the thermostatic state is stable. If lim^^o {X) / (Xq) ~ f{B), the anomalous dimensions, 
caused by quantum turbulence, will affect this operator. If lim^^o {X) / (Xq) ~ <5“" or 
~ exp(a(5“^) for some a, the theory is “renormalizeable” for some observables, while if 
lim^^o {X) / (Xq) ~ or ~ exp(a(5“^) for as that are (X)-specific (one a for the scalar 
and another for the tensor, defined below) the theory is “trivial”, in that taking 5 —)■ 0 
makes the vacuum diverge. In the latter case a limit r^/s —)■ 0 is indeed inconsistent, and the 
degree of divergence could be used to understand the behavior of the limiting rj/s on the 
EoS. The statistics required for a systematic study of important (X) is at the moment well 
beyond our capabilities, but in principle this study is achievable with today’s computing 
technology. 

Scaling in the number of dimensions is also interesting: In 3D, vortices can point in 
random directions, and tend to form an instantaneously disordered “glass” which, over a 
long period of time, seeds the well-known Kolmogorov cascade In 2D, vortices can 

only point in two directions. Given the generally attractive potential between vortices, the 
turbulent phase in 2D is characterized by a 2D crystal of more or less “regular” structures 


33l |. which over long periods of time quenches smaller structures into larger ones, thus 


motivating an “inverse cascade” picture (and effectively negative transport coefficents 34|). 
In ID, vortices are absent altogether, with non-linear corrections being given exclusively by 

ID turbulence 


If the phase we 


sound-waves. As sound waves, unlike vortices, undergo normal scattering 
could be described by the kinetic approaches explored thoroughly in [31 
found can be characterized as fiuctuations-induced turbulence, dimensional scaling will have 
the properties outlined here. 

Given a larger lattice, the structure of the solitonic modes can also be analyzed in space, 
to extract the Prandl width and see how the two phases, if this is indeed what they are, 
vary in observables such as Eqs. fflTl) and ([15]). Eventually, given the rich pure gauge lattice 
QGD data in this area, this opens the door to quantitative exploration of the analogies 
conjectured in [23 1. 

If our picture of a non-trivial “turbulent” vacuum is confirmed, it is still completely 
unclear to what extent this applies for rj/s ^ 0. For a finite viscosity, the partition function of 
Eq. ((7|) will aquire imaginary components, representing the dissipation of collective degrees 
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of freedom into microscopic ones. Lattice techniques described in this work cannot be 
applied to this system, and it is not clear to what extent the extended collective excitations 
are fragile against decohering degrees of freedom. Hence, we still do not know in what 
regime the gradient expansion would fail for the reasons discussed here. The dimensional 
analysis arguments in terms of sound waves given in the Introduction could provide a good 
benchmark, but this is not a rigorous proof. Nevertheless, the effort here can be taken as 
the hrst step to dehne an effective held theory of hydrodynamics that includes both coarse- 
graining parameters: the mean free path and the huctuation scale. Perhaps a mixed classical 
and quantum effective theory 3^ could provide more quantitative answers. 


In conclusion, we found intriguing hints of a non-trivial vacuum structure in ideal hy¬ 
drodynamics without chemical potentials. If conhrmed, it could potentially invalidate the 
gradient expansion as a complete effective theory for a range of microscopic parameters of 
the underlying theory. We hope and expect subsequent work will clarify the properties of 
the non-trivial phase and the transition between the two. 
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